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We study the transfer matrix for domain wall fermions to understand the origin and significance 
of oscillatory contributions to hadron correlation functions that arise for M >1. For a free parti- 
cle in one space, one time, and one flavor dimension, the eigenmodes of the one-body operator 
appearing in the transfer matrix are calculated, and the role of the negative eigenmodes arising 
when M > 1 is studied. In the case of three space dimensions, oscillatory behavior for hadron 
correlation functions in QCD is shown to emerge for free fermions when M exceeds 1, and to 
increase with increasing M. Analogous behavior is observed for domain wall fermions on HYP 
smeared MILC lattices, and a procedure is demonstrated for subtracting oscillating terms from 
physical observables. 
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1. Introduction 

Oscillatory terms in the dependence of lattice correlation functions on Euclidean time t arise 
for domain wall fermions with domain wall mass M > 1, so the goal of this work is to explore 
their origin, whether they are benign or malevolent, and removal of their contributions to physical 
observables. 

For domain wall fermions[|l|, |2|], the domain wall mass M generates a bound state in the extra 
"flavor" direction s. For free fermions, the region < M < 2 is free of doublers, and the choice 
of M < 1 produces exponential decay e^"' in the flavor direction. The choice M > 1 produces 
singularities in a and oscillating behavior for the propagator in the flavor direction!^ G ~ (1 — 
f^^-s-s Pqj. fej-jtiions in a gauge field, an additive mass shift modifies the values M = 0, 1, and 2 
above, and in practice, one uses spectral flow to determine the region of M with the lowest density 
of low eigenmodes of the Hermitian Dirac operator, yjDiy. Thus, for the DBW2 action||3|], the 
LHPC hybrid action using DW fermions on HYP smeared MILC configurations [Q], and the 2+1 
dimensional DW action |^], the domain wall masses were determined to be M = 1.7, 1.7, and 1.8 
respectively, and in the first two cases, oscillatory terms in the t dependence of hadron correlations 
functions have been observed. 

To understand the origin of the oscillatory modes, we study the behavior of the transfer matrix 
int. If one analyzes the lattice theory in the physical dimensions of space and t by integrating out 
all the flavor slices except the boundaries containing the chiral fermions, the theory is intrinsically 
nonlocal and we are unable to construct a one-step transfer matrix. However, we may alternatively 
consider correlation functions in the full dimensions of space, s and t, and restrict our attention 
to operators that are located on the flavor boundaries containing chiral fermions. In this case, it 
is straightforward to construct the DW transfer matix in the ^-direction following the approach 
of Liischerf^, and the analysis is similar to the construction of the DW transfer matrix in the 
s direction[^ |8|]. For simplicity, we will first consider the case of free fermions in one spatial 
dimension plus s plus t and then examine lattice calculations of correlation functions in three spatial 
dimensions plus s plus t for both free fermions and fermions propagating in HYP smeared MILC 
configurations for 2+1 flavors of dynamical fermions. We will show that free fermions and fermions 
in dynamical gluon configurations display the same oscillatory behavior, and implement a simple 
procedure for subtracting the oscillatory contributions from correlation functions. 

2. Domain wall transfer matrix in t 

We write the temporal gauge domain wall action in one spatial dimension in the following 
form: 

DDiy(f,f,5;f',x',5') 

D{xs}{x,s;x ,s) 



X ,s 



B C \ 
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where 



Q(v,v') 



ne{l,s} 

ne{l,s} 

—m, n = s, s = ov Ls — I, 
1, otherwise. 



v=(x,j), z=(?,x), k: = (1,-0, £> = 2. 



The action has the same form in three spatial dimensions, with replaced by Pauh matrices for 
H =x,y,£ and l2x2 for n=s,B multipUed by l2x2, and D = 4. In the free case, the eigenfunctions 
are plane waves and B has the value 



D-l 



B = D + l —M— ^ cos^i — cos^i 
(=1 ^ ' 



<i 



<D-1 



Thus, if M < 1, B is nonnegative, whereas if M > 1, there are values of k such that B is negative. 

A simple way to see the correspondence between the transfer matrix and the standard form of 
the domain wall action is to write the coherent state path integral arising from a one-body Hamilto- 
nian and compare it with the domain wall path integral. We use the following properties of fermion 
coherent states |^): 

Breaking evolution in EucUdean time into steps of size £ in the usual way and combining the terms 
indicated by the braces, we obtain 

t J ^ 

so that the general term in the action connecting times t and ? — 1 is 

S= ^ta^ta + ^ta [^^^"] afi ^(t-l)l}- 

Suppressing details for the lower components for brevity, for the upper components, the general 
term in the lattice action connecting times t and t — I can be brought into this form by writing 

-Sf = Vm^apVtp + Vm^ta{t-l)pV(t-l)P 

= 1/^*6^251/2^^ ^*gl/2g-l/2^g-l/2gl/2^ 

= l/A*Sl/2gl/2^^^*gl/2g^l/2^^^g^gl/2^^ 



3 



Oscillatory terms in the domain wall transfer matrix 



John W. Negele 



B ^-MB 



transforming Grassman variables t, =B^\\f, and identifying e^^^ 

Explicitly constructing the transfer matrix in this way (with appropriate treatment of antiparti- 
cle fields as well) for the domain wall action with one spatial dimension yields 



-H 



1 



B 



-1 



1 Co 
C\ l+c[Co 





Bo 




1 



Tr e 



Hence, to isolate the oscillatory terms in the transfer matrix, we calculate the eigenmodes 

e^^^|f{x,s) = r\^{x,s). (2.1) 



The eigenfunctions in Eq. lA are orthogonal with the weight matrix ( q b) and all the eigenvalues 
are necessarily real. Negative eigenvalues Tj correspond to terms that change sign under evolution 
by one time step and thus produce oscillatory behavior in temporal correlation functions. 



3. Negative eigenmodes for free fermions in 3 dimensions 



The eigenvalues and eigenfunctions in Eq. 2. 1 were computed for free fermions in one space 



plus time and flavor dimensions on a mesh with A''^ = 100, Ns = 20 and ma = 0.05, and typical results 
are shown in Fig. [l[ The eigenfunctions ai^e momentum eigenstates and the quantity E = log\r]\ 
is plotted as a function of p for both positive eigenvalues tj (denoted by red crosses), where it 
corresponds to the energy of the mode, and negative eigenvalues (denoted by blue circles). At M 
= 1, there are only positive modes, as expected, and there is a single physical mode corresponding 
to the free quark dispersion relation. As M increases through 1.1, 1.3, 1.5, and 1.7, the blue points 
show the increasing presence of negative eigenmodes. However, all these modes have energies 
comparable to the positive modes at or above the cutoff scale, and are thus expected to simply be 

additional cutoff effects. The lower right plot shows the chirality R = '^■r'^^+'^j.'^^ ^ which should 

¥r¥r+¥lVl 

be small for chiral modes, for each of the points shown in the middle right plot for M = 1.7. As 
expected, only the free mode is chiral, and the negative modes at the cutoff scale behave the same 
as the positive modes at that scale. Thus all the evidence suggests that these modes are iimocuous 
and can be treated like other cutoff effects. 



4. Oscillations in correlation functions in 5 dimensions 

Analogous negative modes arise for free fermions in three space plus time and flavor dimen- 
sions, and their effect on the proton two-point correlation function C2{t) is shown in the left panel 
of Fig. ^ for a 20^ x 32 lattice with ma = 0.0810 . To display the oscillations, we calculate the 
effective mass Ef,ff = logC{t) /C{t + \) and plot the difference between the effective masses at M 
and M = 1, Eeff{M) —Eeff{\). For clarity, the results at each t are displaced slightly to the right 
with increasing M to avoid overlapping symbols and to indicate the behavior with increasing M. 
Note that oscillations first become visible at M = 1.4 and increase strongly as M is increased to 
1.7, rising on odd time slices and falling on even time slices. 
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0.5 1 1.5 2 



IpJ = 2lsin(:ik/L^)l 

Figure 1: E = log\r\\ as a function of p for positive and negative eigenvalues Tj for free fermions in one 
spatial dimension. The left plots show M = 1, 1.1, and 1.3 and the top right plots show M = 1.5 and 1.7. The 
lower right plot shows the chirality R for all the modes at M = 1.7. 

The same qualitative behavior is observed in the right panel of Fig. § for domain wall fermions 
on HYP-smeared MILC lattices, for which the chopped lattice size is 20^ x 32, a = 0.125 fm, = 
759 MeV and 50 configurations were used. For the MILC lattices, the oscillations become visible 
at M = 1 .5 and grow much larger by M = 1 .7. 

An alternative way of detecting oscillations, which is convenient computationally, is to calcu- 
late the ratio of the correlation function at M to that at M = \, R = C2{M,t)/C2{l,t). This ratio is 
shown in Fig. ^ for the pion and the proton, where for aesthetic purposes, we have also multiplied 
by e^'"' to make the curve flat. Note that oscillations become visible in both the pion and proton 
at M = 1.5, and become large by M = 1.7. Also note that the effect is much larger for the proton, 
where for M = 1.7 it is 12% whereas for the pion it is 3%. 
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Proton mass difference : free 




Proton mass difference : bootstrap over 50 efgs 
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Figure 2: Proton effective mass E^ff = logC{t)/C{t + 1 ) in three space dimensions for free fermions (left) 
and HYP smeared MILC lattices (right) for domain wall masses 1 <M< 1.7. To display the oscillations, 
the plot shows Eeff{M) ~Eeff{l). 
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Figure 3: Correlation function ratio R — CjiM ,t) /C2{1 ,t)e^"" for the pion (left) and nucleon (right) for 
HYP smeared MILC lattices 



5. Subtraction of oscillating terms 

Since all the evidence suggests that the negative eigenmodes are innocuous excitations at the 
cutoff scale for both free quarks and QCD, a straightforward way to remove their effect from 
correlation functions is to fit the correlation function with the sum of an exponential for the physical 
ground state, exponentials for physical excited states, and exponentials multiplied by (—1)' for the 
negative eigenmodes of the transfer matrix. Thus, for the 2-point function, including the leading 
corrections yields 

f{t)=Aoe-'"''+Aie-'"''+Aie-''''i-\y, 
and for the normalized 3-point function for a matrix element we obtain 



(7(|)^(0/(-f)) 



: Bo + cosh[(mi — mo)f] +Bicosh[(mi — mo)f](— 1)'. 



Figure ^ shows the application of this analysis[Q] for MILC lattices with mj^ = 759 MeV. The 2- 
point function determines all three masses as well as three amplitudes, and using these masses, 
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Figure 4: Six-parameter fit to the nucleon two-point function (left) and three-parameter fit to the normalized 
vector current three-point function (right) for DW fermions on HYP-smeared MILC lattices at M = 1.7. 

the normalized 3-point function determines the three amlitudes Bq,Bi, and Bi and thereby enables 
accurate extraction of the physical matrix element, Bq, for the vector current operator. 

6. Conclusions 

In summary, we have calculated the transfer matrix in the Euclidean time direction for free 
domain wall fermions and shown that negative eigenmodes arise for M > 1 and produce oscillatory 
behavior in lattice correlation functions. We argue that these are innocuous lattice effects at the 
cutoff scale and show how their effects can be removed from QCD correlation functions. In the 
future, however, since the additive mass shift depends on the gauge action and lattice spacing, we 
believe care should be taken in choosing the domain wall mass M for each action and lattice spacing 
to eliminate the appearance of these modes in lattice calculations. 

This work was supported in part by funds provided by the U.S. Department of Energy un- 
der grant DE-FG02-94ER40818 and computer resources were provided on the MIT Blue Gene 
computer. 
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